High-pressure phases of silane 
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It has been suggested that hydrogen may metallise at lower pressures if it is "precompressed" . 
Here we introduce a search strategy for predicting high-pressure structures and apply it to silane 
using first-principles electronic structure computations. It is based on relaxing randomly chosen 
structures, and is demonstrated to work well for unit cells containing up to at least ten atoms. 
We predict that silane will metallise at higher pressures than previously anticipated, but we suggest 
that the metallic phase might show high-temperature superconductivity at experimentally accessible 
pressures. 

PACS numbers: 61.66.Fn,71.20.-b,74.10.-|-v 
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As early as 1935 Wigner and Huntingdon suggested 
that hydrogen would become metallic under sufficient 
compression. |lj There have been many attempts to pro- 
duce metallic hydrogen under high-pressure laboratory 
conditions, but it has stubbornly remained insulating in 
diamond-anvil-cell experiments up to the highest pres- 
sure so far achieved of 342 GPa.|3| It has recently been 
suggested|^ that the group IV hydrides, methane (CII4), 
silane (SiH4), and germane (GeH4), might become metal- 
lic at pressures achievable in diamond anvil cells because 
the hydrogen in these materials is "chemically precom- 
pressed" by the presence of the group IV atoms. Such 
metallic phases might exhibit phonon-mediated high- 
temperature superconductivity. 

Recently, Feng et aZ.|j| reported a first-principles 
density-functional-theory (DFT) study of silane at high 
pressures. They predicted that a layered structure would 
be stable at pressures above 25 GPa. This phase was 
found to be semi-metallic at 91 GPa, and the density of 
electronic states at the Fermi energy increased rapidly 
with further application of pressure. An estimate us- 
ing the Bardcen-Cooper-Schrieffer (BCS) theory [5| in- 
dicated that the superconducting transition temperature 
would increase from about K at 91 GPa to about 166 K 
at 202 GPa, pressures which are comfortably within 
range of diamond anvil cells. 

Accurate cohesive energies and equilibrium volumes 
for materials of known crystal structure have been ob- 
tained using first-principles methods, but the prediction 
of stable structures has remained very difficult. Using 
these methods, it is routine to relax a structure from a 
chosen initial configuration to a nearby minimum in the 
potential energy surface. Relaxing from a set of initial 
configurations may lead to several distinct local-energy 
minima, the lowest-energy one corresponding to the most 
favoured structure. However, if none of the initial con- 
figurations is sufficiently close to the global minimum 
then the configurations will become trapped in higher- 
energy local minima. Trapping in local-energy minima 
is a severe problem because the number of minima is 



expected to increase exponentially with the number of 
degrees of freedom. In the past this problem has typ- 
ically been addressed by selecting initial configurations 
derived from experimental information known about the 
system in question, the known structures of similar ma- 
terials, chemical intuition, and results from simpler com- 
putational methods. The shortcomings in this approach 
have led to a growth in using first-principles methods 
in combination with more adva_nced_search techniques to 
determine stable structures. 



We have developed a simple strategy for generating ini- 
tial configurations, which consists of choosing essentially 
random unit cell parameters and atomic positions. We 
first choose the number and identity of the atoms in the 
unit cell. To form "random" unit cells we choose three 
cell-vector lengths randomly and uniformly between 0.5 
and 1.5 (in arbitrary units), and three cell angles ran- 
domly and uniformly between 40° and 140°. The result- 
ing cell vectors are scaled to produce the desired volume, 
which is also selected randomly and uniformly from be- 
tween 0.5 and 1.5 of some chosen volume. The positions 
of the atoms are then chosen randomly and uniformly 
within the cell. First-principles methods are then used 
to relax each structure until a minimum in the enthalpy 
is reached. The structures obtained from such a proce- 
dure are clearly liable to be trapped in local minima, and 
undoubtedly the probability of finding the global mini- 
mum will drop rapidly with increasing complexity of the 
system. However, for systems of moderate complexity, 
tests indicate that our approach has a high probability 
of locating the global minimum. In practice, the number 
of initial configurations required to find the global min- 
imum depends on the system, and particularly on the 
number of degrees of freedom it possesses. Our approach 
has been to continue generating configurations until the 
relaxed structures with low energies are generated several 
times and, where possible, to look for the occurrence of 
previously-known "marker" structures. 

All of the first-principles calculations were performed 
using a developer's version of the CASTEP code.yi| For 



the initial search over structures we used the local den- 
sity approximation (LDA) for the exchange-correlation 
functional [Il|, and default ultrasoft DseudoDotentials.[l^ 
The fc-point sets were generated separately for each unit 
cell encountered during the procedure, and no symme- 
try restrictions were applied at any point. Medium 
quality Brillouin zone sampling using a grid of spac- 
ing 27rx 0.07 A~^ and a plane wave basis set cutoff of 
280 eV were found to be sufficient for the initial search 
over structures. When re-calculating the enthalpy curves 
with higher accuracy we used a Brillouin zone sampling of 
27rx 0.03 A~^, harder pseudopotentials requiring a plane 
wave cutoff of 360 eV, and the Perdew-Burke-Ernzerhof 
(PBE) Generalised Gradient Approximation (GGA) den- 
sity functional 13] to aid comparison with Ref. [j|. Each 
of the calculations was for the zero-temperature ground 
state and we neglected the zero-point vibrational energy. 

We have applied our approach to silicon at zero pres- 
sure. Using two-atom unit cells and 100 randomly chosen 
initial configurations, we obtained the global minimum 
energy diamond structure (15 times), f3-tm (2 times), 
Imma (14 times), and simple hexagonal (15 times). Both 
experiment and DFT calculations agree that these phases 
are the four lowest-pressure thermodynamically stable 
phases of silicon. jlj| A study of 8-atom unit cells of car- 
bon at 300 GPa generated the diamond structure as the 
lowest-enthalpy phase in addition to numerous higher- 
enthalpy phases, including the BC8 phase, which has 
been postulated as a stable phase of carbon at TPa 
pressures, [l^ We also studied 16-atom unit cells of hy- 
drogen at 250 GPa, finding, among other phases, the 
Cmca molecular phase which DFT calculations have pre- 
dicted to be the most stable at that pressure. p?7j 

Our approach is particularly suitable for the predic- 
tion of the structures of high-pressure phases because 
they often have fairly simple structures with a small 
number of atoms in the primitive unit cell. In addi- 
tion, our method is straightforward and does not re- 
quire the selection of highly-system-specific parameter 
values. Our scheme could be modified to include other 
features, such as methods for extricating configurations 
from local energy minima. |l8l| However, such modifica- 
tions would complicate the method and add to the com- 
putational cost, which would limit the number of initial 
configurations which can be sampled. It is also possible 
to introduce other constraints on the initial configura- 
tions so that one can predict the structures of clusters, 
defects, surfaces, etc. For example, we have used our 
approach to search for low-energy self-interstitial defect 
structures in silicon. To construct the initial configura- 
tions we took a 32-atom body-centred-cubic unit cell of 
diamond-structure silicon, removed an atom and its four 
nearest neighbours, and then placed six silicon atoms ran- 
domly within a spherical region of radius about 1.5 bond 
lengths centred on the first atom that we removed. Us- 
ing our approach, we readily found the split-(llO) and 



hexagonal interstitial configurations, which other DFT 
calculations have shown to be lowest in energy [l9|, and 
also predicted that all other self-interstitial local energy 
minima are at least 0.3 eV higher in energy. The genera- 
tion of numerous low-energy structures is an interesting 
and useful feature of our strategy, as they can give in- 
formation about the structures which might be formed 
under other conditions of pressure and temperature, or 
out of thermal equilibrium. 

We tested our search strategy for simulation cells con- 
taining one, two, and four SiH4 units at a pressure of 
250 GPa. With one SiH4 unit per cell the lowest-enthalpy 
structure occurred four times from 13 initial configura- 
tions, but it was 0.36 eV per SiH4 unit higher in enthalpy 
than the most stable structure found with two SiH4 units 
per cell, which occurred 8 times from 37 initial configu- 
rations. With four SiH4 units per cell the most stable 
structure was 0.13 eV per SiH4 unit higher in enthalpy 
than for the two SiH4 unit cell, and each of the low- 
enthalpy phases occurred only once from 38 configura- 
tions. Clearly the simulation cells containing one SiH4 
unit are too constrained to produce very-low-enthalpy 
structures, while those containing four Sill4 units have 
too many degrees of freedom for our search strategy to 
find the global minimum enthalpy structure with the 
number of initial configurations used. We concluded that 
using two Sill4 units per simulation cell and 40 initial 
configurations represented an affordable and efficient ap- 
proach to finding low-enthalpy phases of silanc at high 
pressures. 

Further runs were performed at 0, 50, 100, 150, and 
200 GPa, with two Sill4 units per simulation cell and 
40 initial structures. At GPa the lowest-energy struc- 
ture could not be determined with confidence because the 
low-energy structures were not found repeatedly. This re- 
flects the fact that at zero pressure many different struc- 
tures have similar energies (within 0.1 eV per SiH4 unit), 
including those consisting of silane molecules, disilane 
plus a hydrogen molecule, and various polymeric forms 
plus hydrogen molecules. The various packings of silane 
molecules were therefore inadequately explored using just 
40 initial structures. The lowest-energy structure found 
(by about 40 meV per Sill4 unit) was a low-symmetry 
silane molecular crystal which was about 50 meV per 
SiH4 unit less stable than the higher-symmetry Tl pack- 
ing considered by Feng et al. |j|. 

Having selected the low-enthalpy phases from the 
above runs, we re-calculated their enthalpies with greater 
computational accuracy. The calculated enthalpies of 
various structures are shown as a function of pressure in 
Fig. n The details of some of the lower-enthalpy struc- 
tures found are described in the auxiliary material to this 
Letter. |l5!| 

An insulating chain-like structure with P2/c symme- 
try was found to be the most stable in a small region of 
pressures around 40 GPa, but a structure of /4i/a sym- 
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FIG. 1: The enthalpies per SiH4 unit of various structures as 
a function of pressure, referenced to the Tl phase of Ref. 3 . 
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FIG. 2: The densities of states of the /4i/a and C2/c phases 
at their predicted coexistence pressure of 262.5 GPa. The 
partial density of states projected onto the H atoms is shown 
by dashed curves, and the Fermi energies are indicated by the 
vertical dashed lines. 



metry became more stable at about 50 GPa. The lii/a 
phase is predicted to be insulating at 50 GPa but, as the 
pressure is increased, the bandgap gradually closes and 
it becomes semi- metallic, as can be seen in the density of 
states plotted in Fig. |2| Given the well-known shortcom- 
ings of standard DFT approaches such as the LDA and 
GGA functionals, which normally underestimate band 
gaps, it is likely that our calculations underestimate the 
pressure at which the /4i/a phase becomes semi-metallic. 
The /4i/a phase remains the lowest-enthalpy phase up to 
262.5 GPa, at which pressure a denser structure of C2/c 
symmetry becomes stable. This phase is a good metal, 
see Fig. El 

In the /4i/a structure, which is illustrated in Fig. 13 
each Si atom is bonded to 8 H atoms. Each of the 
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FIG. 3: The /4i/a and C2/c structures. The golden spheres 
represent silicon atoms, and the white spheres hydrogen. Note 
the Si2H2 planes in the Mi /a structure, and the two, three 
and four-fold coordinated H atoms in the C2/c structure. 



H atoms forms a "bridge" between two neighbouring 
Si atoms. There are two such bridges between every 
pair of neighbouring Si atoms, the four atoms forming 
a Si2H2 plane. This bonding arrangement consists of 
two "banana bonds" , which are electron-deficient three- 
centre-two-electron bonds, reminiscent of those in B2H6 
(diborane).[20| The arrangement of the Si2H2 planes 
keeps the H atoms away from each other, and the /4i/a 
structure contains only Si-H bonds. Each of the Si sites 
and the H sites is equivalent in this high-symmetry struc- 
ture which has a 10-atom primitive unit cell. The lay- 
ered 03 structure proposed in Ref. 4] also contains some 
bridging Si-H-Si bonds, although it also contains H atoms 
which are bonded to single Si atoms. Our /4i/a struc- 
ture is about 0.4 eV per SiH4 unit lower in enthalpy than 
the 03 structure at 100 GPa, and about 0.7 eV per SiH4 
unit lower at 180 GPa. 

In the metallic C2/c structure illustrated in Fig. O 
each Si atom is bonded to 11 H atoms. Each of the 
Si sites is equivalent and there are three inequivalent H 
sites in which the H atom forms bonds to two, three, or 
four Si atoms. The average coordination number of the 
H atoms in the C2/c structure is larger than in /4i/a, 
which accounts for its smaller equilibrium volume. The 
metallic character of this phase means that the nature 
of the inter-atomic bonding is less clear cut than in the 
14:1 /a structure, but again a description in terms of elec- 
tron deficient Si-H bonds is appropriate. 

Our search also revealed a number of other low- 
enthalpy phases, including the 03 Pmna structure stud- 
ied by Feng et aLjJ] We found that, at around 200 GPa, 
03 Pmna converts into a more compressed structure of 
the same symmetry, but it still has a substantially higher 



enthalpy than our Mi /a phase. We followed the more 
compressed phase to lower pressures and found it to be 
more stable than the 03 Pmna structure down to about 
110 GPa. We also found a low-enthalpy phase with IA2d 
symmetry which is only 0.1 eV per SiH4 unit above the 
most stable phase at 100 GPa. The P2/c phase men- 
tioned above converts to a structure with PA/nbm sym- 
metry at 100 GPa. Interestingly, the PA/nbm phase has 
a layered structure, like the 03 structure, but with every 
hydrogen atom participating in a bridging bond, which 
results in its lower enthalpy. 

Having found the low-enthalpy metallic C2/c structure 
using simulation cells containing two SiH4 units, we stud- 
ied its stability to rearrangements of the H atoms using 
simulation cells containing four SiH4 units at 250 GPa. 
To make initial configurations we used the C2/c Si frame- 
work and placed the 16 H atoms randomly within the 
simulation cell, and then relaxed all of the atomic posi- 
tions and cell parameters. Using 15 such initial config- 
urations we obtained the C2/c structure as the lowest- 
enthalpy phase, an Fddd structure only 30 meV per SiH4 
unit higher in enthalpy, and several other slightly less 
stable phases. All of the low-enthalpy phases contained 
a mixture of two-, three-, and four- fold coordinated H 
atoms, as in the C2/c structure, and they all had very 
similar volumes and densities of electronic states. 

Silane is in fact unstable to decomposition into hydro- 
gen and silicon at zero pressure, and we therefore con- 
sidered the thermodynamic stability of our new silane 
phases against processes involving the formation of hy- 
drogen molecules. At low pressures we found some low- 
enthalpy structures containing covalently bonded hydro- 
gen molecules, indicating the instability of silane at these 
pressures, but none of the low-enthalpy structures at high 
pressures contain any H-H bonds, suggesting stability 
against decomposition. We investigated this further by 
considering decomposition into diamond-structure Si and 
molecular hydrogen at zero pressure, into simple hexag- 
onal silicon [lj| and Cmca-structure hydrogen [l2| phases 
up to 40 GPa, and into face-centred cubic silicon and 
Cmca-structure hydrogen at higher pressures. While the 
03 structure proposed in Ref. Q is found to be unsta- 
ble with respect to decomposition into face-centred cu- 
bic silicon and Cmco-structure hydrogen at its metalliza- 
tion pressure, our /4i/a and C2/c structures are stable 
with respect to such a decomposition. This lends fur- 
ther weight to our contention that the C2/c structure of 
silane, or something very like it, is thermodynamically 
stable at high pressures. 

We now discuss the possibility of phonon-mediated su- 
perconductivity in the C2/c phase. Within the BCS 
theory, the superconducting transition temperature de- 
pends on the values of three parameters: the density of 
states at the Fermi energy of the normal state, the Debye 
temperature of the phonons, and the effective electron- 
phonon coupling parameter. We estimated the densities 



of states at the Fermi energies of the C2/c and /4i/a 
phases at their coexistence pressure of 262.5 GPa using 
DFT methods. As shown in Fig. |21 the C2/c phase is 
a good metal, with a density of states at the Fermi en- 
ergy of 1.6 X 10~^ eV~^ A, which is about 46% of the 
free-electron value and is close to the zero-pressure value 
for lead. We decomposed the densities of states into H- 
and Si-associated components using a MuUiken popula- 
tion analysis. [2]| Fig. El shows that below the Fermi en- 
ergy the density of states is a little larger on the H atoms 
than on the Si ones, while above the Fermi energy it be- 
comes progressively larger on the Si atoms. This is con- 
sistent with the fact that the Pauling electronegativity of 
H (2.2) is a httle larger than that of Si (1.9). The fact 
that the density of electronic states at the Fermi energy 
projected onto the H atoms is large suggests that the cou- 
pling between the electrons close to the Fermi energy and 
phonon modes involving motions of the H atoms may be 
significant. We also performed DFT calculations of the F- 
point phonon frequencies of the C2/c phase at 262.5 GPa, 
finding the highest phonon frequency to be 2500 cm~^. 
This corresponds to an estimated Debye temperature of 
3600 K, which is about 40 times larger than that of lead 
at zero-pressure. An examination of the phonon eigen- 
vectors indicates that the highest-frequency modes cor- 
respond largely to motions of the H atoms. The very 
large Debye temperature and the substantial density of 
states at the Fermi energy suggests that the C2/c phase 
might exhibit high-temperature superconductivity, but 
one should be cautious about this conclusion. A reliable 
estimate of the superconducting transition temperature 
within the BCS theory also requires evaluation of the 
full electron-phonon coupling matrix elements, which is 
beyond the scope of this work. 

In summary, we have used a random searching strategy 
in conjunction with first-principles electronic structure 
computations to predict the stable high-pressure phases 
of silane. We find insulating/semi-metallic behaviour up 
to 260 GPa, at which point a first-order phase transi- 
tion occurs to a good metal. This phase might show 
high-temperature superconductivity at a pressure which 
is achievable within a diamond anvil cell. 
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phonon calculations. 
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of accuracy described in the letter. Only the fractional coordinates of symmetry inequivalent atoms are reported. 



017006 (2006). [13] 

[5] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. 

Rev. 108, 1175 (1957). [14] 

[6] R. Martonak, A. Laio, and M. Parrinello, Phys. Rev. 

Lett. 90, 075503 (2003). [15] 

[7] S. Yoo and X. C. Zeng, Angewandte Chemie Interna- [16] 
tional Edition 44, 1491 (2005). [17] 

[8] S. Gocdccker, W. Hellmann, and T. Lcnosky, Phys. Rev. 

Lett. 95, 055501 (2005). [18] 

[9] A. R. Oganov, C. W. Glass, and S. Ono, Earth and Plan- 
etary Science Letters 241, 95 (2006). [19] 
[10] M. D. Segall, P. J. D. Lindan, M. J. Probert, C. J. [20] 
Pickard, P. J. Hasnip, S. J. Clark, and M. C. Payne, 
J. Phys: Condens. Matter 14, 2717 (2002). [21] 
[11] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 

(1981). 
[12] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 



J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 

A. Mujica, A. Rubio, A. Munoz, and R. J. Needs, Rev. 

Mod. Phys. 75, 863 (2003). 

See EPAPS document No. xxxxxxxxxxxx. 

M. T. Yin, Phys. Rev. B 30, 1773 (1984). 

K. A. Johnson and N. W. Ashcroft, Nature 403, 632 

(2000). 

D. J. Wales, Energy Landscapes (Cambridge University 

Press, 2003), chap. 6. 

R. J. Needs, J. Phys: Condens. Matter 11, 10437 (1999). 

L. Burnelle and J. J. Kaufmann, J. Chem. Phys. 43, 3540 

(1965). 

M. D. Segall, R. Shah, C. J. Pickard, and M. C. Payne, 

Phys. Rev. B 54, 16317 (1996). 



